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Abstract 

I review recent and not so recent progress on formulating and numerically im¬ 
plementing a consistent set of relativistic equations which describe the space- 
time evolution of viscous relativistic fluids without violating causality. 


1 Introduction 

Ideal fluid dynamics has been used successfully to predict the collective flow patterns in Au+Au colli¬ 


sions at the Relativistic Heavy-Ion Collider RHIC (for a review see The ideal fluid description works 


well in almost central Au+Au collisions near midrapidity at top RHIC energy, but gradually breaks down 


in more peripheral collisions, at forward rapidity, or at lower collision energies (2}. To describe such 


deviations from ideal fluid dynamics quantitatively and use the experimental data to extract values or 
phenomenological limits for the transport coefficients of the hot and dense matter created during the 
collision requires the numerical implementation of dissipative relativistic fluid dynamics. Although a 
formulation of such a theory which avoids the longstanding problems of acausal signal propagation and 
other instabilities associated with the original relativistic fluid equations given by Eckart and Landau & 
Lifshitz, has been known for almost 30 years o, significant progress towards its numerical implemen¬ 
tation has only been made very recently Ena mm. At this point, we are only at the very beginning of a 
program that will eventually apply a viscous relativistic fluid dynamical approach to heavy-ion collision 
data. Existing numerical implementations are (l+l)-dimensional and can only describe cylindrically 
symmetric transverse expansion with boost-invariant longitudinal dynamics (6j[7). The codes are still in 
the process of being tested. In this contribution I try to give a pedagogical summary ® of the theoretical 
background and summarize a few first results. 

2 Ideal fluid dynamics 

Before explaining the structure of the equations for causal dissipative relativistic fluid dynamics, let me 
quickly review the case of ideal fluid dynamics. Any fluid dynamical approach starts from the conserva¬ 
tion laws for the conserved charges and for energy-momentum. 


d,jN^ = 0, i = l,...,k, 

d ljT ' JV = 0. 


( 1 ) 

( 2 ) 


For simplicity I will restrict myself to k = 1 (say, N lt = net baryon number current) and drop the index i 
in 0 . It must also ensure the second law of thermodynamics 


d^ > 0 


( 3 ) 


where S 1 '' is the entropy current. Ideal fluid dynamics follows from these equations under the assumption 
of local thermal equilibrium, i.e. if the microscopic collision time scale is very much shorter than any 


macroscopic evolution time scale such that the the underlying phase-space distribution f(x,p ) relaxes 


essentially instantaneously to a local equilibrium form (upper signs for fermions, lower signs for bosons) 
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( 4 ) 
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where vJ' ix) is the local fluid velocity at point x, p(x) is the local chemical potential associated with the 
conserved charge N (it enters with opposite sign in the distribution / for antiparticles), and T(x) is the 


local temperature. Plugging this into the kinetic theory definitions 

N^(x) = J ^-p' i [f(x,p) - f(x,p)], (5) 

T^(x) = J < ^-p tJ, p v [f(x,p) + f(x,p)\, (6) 

S^( x ) = - j ^p^if (x,p) In f(x,p) ± (l=F/(x,p))ln(l=F/(x,p)) + (/ <-► /)] (7) 

leads to the ideal fluid decompositions 

= nu>\ ( 8 ) 

T% = e//-pA^ (with A ^ = g^-u»u u ), (9) 

S eq = SU ^i (! 0 ) 


where the local net charge density n, energy density e, pressure p and entropy density s are given by the 
standard integrals over the thermal equilibrium distribution function in the local fluid rest frame and are 
related by the fundamental thermodynamic relation 


T s = p — pn + e. (11) 

Inserting Eqs. <15171 into Eqs. CD yields the relativistic ideal fluid equations shown in Eqs. can 
below. Using Eq. CD together with the Gibbs-Duhem relation dp = s dT + n dp, it is easy to prove 
that, in the absence of shock discontinuities, these equations also conserve entropy, i.e. 8^ S 1 ' = 0. 

Note that the validity of the decompositions <151171 only requires local momentum isotropy (i.e. 
that in the local fluid rest frame the phase-space distribution reduces to a function of energy E only, 
f(x,p) = f (p ■ u (x )'.T{ x), /i (x))), but not that the distribution function has the specific exponential 
form 0 which maximizes entropy. This may have relevance in situations where the time scale for local 
momentum isotropization is much shorter than for thermalization (i.e. it is much easier to change the 
direction of the particles' momenta than their energies), with the macroscopic hydrodynamic time scale 
in between. In this case the local microscopic states would not maximize entropy, and the relation (CD 
would not hold between the quantities e, p, n, and s defined through eqs. (Blflol . Still, they would 
follow ideal fluid dynamical evolution since entropy production by microscopic kinetic energy-shifting 
processes would happen only on time scales which are large compared to the macroscopic evolution time 
scales. 

The ideal fluid equations read (with 9 = 8-u denoting the local expansion rate) 

n = 
e = 


where we decomposed the partial derivative d' 1, 
components D = u v d v and V M = t\ >w d v which in the local fluid rest frame reduce to the time derivative 
f = Df and spatial gradient. The first two equations describe the dilution of the local baryon and energy 
density due to the local expansion rate 9, while the third describes the acceleration of the fluid by the 
spatial (in the local frame) pressure gradients, with the enthalpy e+p acting as inertia. The 5 equations 
(TT2IT41) for the 6 unknown functions n, e, p, u 1 ' (remember that vJ l u ll: = 1) must be closed by supplying 
an Equation of State (EOS) p = p(n, e). 


-(e + p) 9, 

\/Up 

e + p ’ 


(13) 

(14) 


= u 1 ' I) + V 1 ' into “longitudinal” and “transverse’ 
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3 Non-ideal fluid decomposition 

As the hydrodynamic evolution changes the local energy and baryon density, microscopic process at¬ 
tempt to readjust the local phase-space distribution to corresponding new local temperatures and chem¬ 
ical potentials. If this does not happen sufficiently fast, the phase-space distribution will start to deviate 
from its local equilibrium form ©: f(x,p) = f eq (p-u(x); T(x), p(x)) + 6f(x,p). The optimal values 
for the (readjusted) local temperature and chemical potential in the first term are fixed by imposing the 
“Landau matching conditions” 

u tl 5T ,w u v = J (u-p) 2 5f(x,p) = u^5N^ = J ^ u-pSf(x,p ) = 0. (15) 


The remaining deviations Sf from local equilibrium generate additional terms in the decompositions of 
N 1 - 1 , T IW , and S'b 

= JV£ + SN^ =n/ + V», (16) 

T^ + 5T^ U = eu^ l u u - {p + U)A^+ tt^+ W fl u u + W u u^, (17) 

S» = S^ + SS^ = nu^ + & 1 . (18) 


The new terms describe a baryon number flow V 1 ' = A lllJ N u in the local rest frame, an energy flow 
W 11 = V ,J ’ + q 1 ' (where q v is the “heat flow vector”) in the local rest frame, the viscous bulk pressure 


II = -iAm,T / "' — p (which contributes to the trace of the energy momentum tensor), the traceless 


T ra (where the 


viscous shear pressure tensor it pu = T^ = 2 (A ll(J A VT +A ,MJ A ,J ' T ) — ^A^A 0 
expression (pis) is a shorthand for “traceless and transverse to u /( and u u ”, as defined by the projector in 
square brackets), and an entropy flow vector ( \> IJ in the local rest frame. 

Inserting the decompositions mm into the conservation laws .rmti yields the non-ideal fluid 
equations 


n 

e 

(e+p+n) vfl 


-nd- V-V + V-u, 

-(e+p+n) 0 + 7- V-W + 2 W-u, 


V^(p+n) - + tt ^ u u 


a^w,, + w^e + (w-v)u^ 


(19) 

( 20 ) 
( 21 ) 


The matching conditions G3 leave the choice of the local rest frame velocity ii 1 ' ambiguous. This 
ambiguity can be used to eliminate either V 1 ' from Eq. (1 1 61 (“Eckart frame”, no baryon flow in the local 
rest frame), in which case the energy flow reduces to the heat flow vector or from Eq. o 

(“Landau frame”, no energy flow in the local rest frame), in which case there is a non-zero baryon flow 
V 11 = — pq-^/1 1 ' due to heat conduction in the local rest frame. (Intermediate frames are also possible, but 
yield no practical advantage.) For systems with vanishing net baryon number (as approximately realized 
in RHIC collisions) the Eckart frame is ill-defined, so we will use the Landau frame. In this frame, for 
baryon-free systems with n = 0 and no heat conduction, the non-ideal fluid equations mm simplify to 

e = — (e+p+n) 9 + n^ u u \ (22) 

(e+p+n) = V^(p+n) - A^V CT 7T W + TT^ v u v . (23) 


The non-equilibrium decompositions (Em involve l+3+5=9 additional dynamical quantities, 
the “dissipative flows” II, q 11 , and tt iiu (the counting reflects their transversality to u 1 ' and the traceless¬ 
ness of 7 t^ v ). This means that we need 9 additional dynamical equations which should be compatible 
with the underlying transport theory for the non-equilibrium deviation Sf(x,p). For the baryon-free case 
without heat conduction, the number of needed additional equations reduces to 6. 
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4 Kinetic equations for the dissipative flows 

The key property of the kinetic equation governing the evolution of the phase-space distribution function 
/ = /eq+4/ is that the collision term satisfies the second law of thermodynamics < 0 , i.e. entropy is 
produced until the system has reached a new state of local thermal equilibrium. We don’t want to solve 
the kinetic theory; instead, we want to write down a phenomenological macroscopic theory which is 
consistent with the constraints arising from the underlying kinetic theory, in particular the 2 nd law. The 
macroscopic theory will be constructed from an expansion of the entropy production rate in terms of 
the dissipative flows which themselves are proportional to the off-equilibrium deviation 8f of the phase- 
space distribution. Assuming the latter to be small, |<5/| <C /,. q |. this expansion will be truncated at some 
low order in the dissipative flows SN 1 ', 8T ,iV . The expansion will involve phenomenological expansion 
coefficients which, in principle, should be matched to the kinetic theory ED In practice, they will often be 
considered as phenomenological parameters to be adjusted to experimental data. In the end, the extracted 
values must then be checked for consistency with the entire approach, by making sure that the dissipative 
corrections are indeed sufficiently small to justify truncation of the expansion a posteriori. 

The equilibrium identity CEB can be rewritten as 

•Seq = P(a, - aN? q + p„Tg , (24) 

where a = ^ • &= p- and (3 U = ^. The most general off-equilibrium generalization of this is 0 

S» = S£ij + = p(a, /3)/3 M - aN^ + (3 V T U ^ + 5T^), (25) 

where, in addition to the first order contributions implicit in the second and third terms of the r.h.s., Q 1 ' 
includes terms which are second and higher order in the dissipative flows SN^ and ()T ,W . [Note that, by 
using the identity CEB between the equilibrium quantities, Eq. CZB can be written in the simpler-looking 
form = s vJ'+^+Q 1 ' but this is not helpful for calculating the entropy production rate.] 

The form of the expansion <EB is constrained by the 2 nd law d^S 1 ' > 0. To evaluate this constraint 
it is useful to rewrite the Gibbs-Duhem relation dp = s dT + n dp as 

dp (p(a, = N^d^a - d^ v . (26) 

With additional help from the conservation laws (1 1 12I) . the entropy production then becomes 

8^ = -5N^d^a + + duQ*. (27) 

Using Eqs. mm to express 5N^ and 6T ,,V in terms of the scalar, vector and tensor dissipative flows 
II, q 1 '. and and introducing corresponding scalar, vector and tensor thermodynamic forces (in 
terms of gradients of the thermodynamic equilibrium variables) which drive these dissipative flows, 
X = -6 = -V-u, X u = ^ - u u = - and X^ = X^u v > (note that X^=X^ is 

traceless and transverse to u), the 2 nd law constraint can be further recast into 

TdpS* = UX - q^X^ + ^ V X^ + Td^Cf > 0. (28) 

Note that the first three terms on the r.h.s. are first order while the last term is higher order in the 
dissipative flows. 


4.1 Standard dissipative fluid dynamics (first order theory) 

The standard approach (which can be found, for example, in the Fluid Dynamics volume of Landau & 
Lifshitz) one neglects the higher order contributions and sets Q /l = 0. The inequality m can than be 
always satisfied by postulating linear relationships between the dissipative flows and the thermodynamic 
forces, 


nT 2 

q u = —X—— V 1 ' 
e+p 


n = -to, 


a 

T 


7= 2 rjV {tl u u) 


(29) 
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with positive transport coefficients Q > 0 (bulk viscosity ), A > 0 (heat conductivity), and rj> 0 (shear 
viscosity ): 


Td-S 


n 2 

T 


q a qa ^^0 > 0 
2 AT 2r/ - ' 


(30) 


(The minus sign in front of the second term is necessary because f, being orthogonal to u 1 '', is spacelike, 
q 2 < 0.) These are the desired 9 equations for the dissipative flows. 

Unfortunately, using these relations in the hydrodynamic equations (USED leads to hydrodynamic 
evolution with acausal signal propagation: if in a given fluid cell at a certain time a thermodynamic force 
happens to vanish, the corresponding dissipative flow also stops instantaneously. This contradicts the 
fact that the flows result from the forces through microscopic scattering which involves relaxation on a 
finite albeit short kinetic time scale. To avoid this type of acausal behaviour one must keep Q 1 '. 


4.2 Second order Israel-Stewart theory 

A causal theory of dissipative relativistic fluid dynamics is obtained by keeping Q 1 ' up to terms which are 
second order in the irreversible flows. For simplicity I here consider only the baryon-free case n = = 0, 

for a general treatment see (20. One writes |0 

cr = - (/3 0 n 2 + /Wtt^) ^ (31) 

(with phenomenological expansion coefficients /3o, Ih) and computes (after some algebra using similar 
techniques as before) the entropy production rate as 


Td-S = n 


Poll - UTdn 



+ 7T *P 

V 2 T )\ 



^(a H0) P2^a0 



(32) 


From the expressions in the square brackets we see that the thermodynamic forces —9 and are 

now self-consistently modified by terms involving the time derivatives (in the local rest frame) of the 
irreversible flows II, it a p. This leads to dynamical (“transport”) equations for the latter. We can ensure 
the 2 nd law of thermodynamics by again writing the entropy production rate in the form (without 
the middle term), which amounts to postulating 


ri 


^a0 


1 

Ai 

1 

Tn 


n + C 0 + ncT^(^0] p 

Ka/3 - 2 yV {a u 0) + Tt^rjTd^ 


i fn + C0], 

T n 


/ r^V 

1 

\ 2 r/T j _ 

Tr 


^a/3 2 (a u /3) 


(33) 

(34) 


Here I replaced the coefficients /3o,2 by the relaxation times r n = fif and t- = 2r//T • In principle both 
r] and t u ,t~ should be calculated from the underlying kinetic theory. We will use them as phenomeno¬ 
logical parameters, noting that for consistency the microscopic relaxation rates should be much larger 
than the local hydrodynamic expansion rate, ^ < 1. 

Let me shortly comment about the approximation in the second equalities in Eqs. d33l34li : We 
are using an expansion scheme for the entropy production rate in which the thermodynamic forces and 
irreversible flows are assumed to be small perturbations. The approximation in Eqs. (1331341 neglects 
terms which are products of the irreversible flows with gradients of the thermodynamic equilibrium 
quantities which are of the same order as the thermodynamic forces. These terms are thus effectively 
of second order in small quantities and should, for consistency, be neglected relative to the other terms 
in the square brackets which are of first order. If one wants to keep them (as done by Muronga 00), 
one should also keep third-order terms in the entropy flow vector Cf‘ for consistency. Of course, where 
the thermodynamic forces and irreversible flows are really small, it shouldn’t matter whether we keep or 
drop these terms. In practice, however, one will use this approach when dissipative effects are expected 
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to be significant, and the dropped terms may not be extremely small. In this case I believe that dropping 
them is more consistent than keeping them. 

There is another reason for dropping these terms: without them, Eqs. d33l34l> are relaxation equa¬ 
tions which describe (in the local rest frame) exponential relaxation (on the time scales 77 n ) of the 
irreversible flows to the values given by Eqs. cn in the first order theory. However, if these terms are 
kept, one has instead equations of the form 


ri = 

[n + (9 + nc 7n ] 

!+7 n C 

n+ c e 

1 

u + e'e 

T n 


T n 

l l+7nC J 

r n 

J 


and similarly for the shear pressure tensor. One sees that both the kinetic relaxation time and the viscosity 

(T \ 

j which involves the macroscopic expansion rate <) ll u 11 . This 
violates the intuition that these transport coefficients should be expressible through integrals of the kinetic 
collision term which involve only microscopic physics (cross sections, local densities, etc.) 

In the second order Israel-Stewart formalism, one thus solves the dissipative hydrodynamic equa¬ 
tions GUEB simultaneously with the kinetic relaxation equations (I33I34I> for the irreversible flows. Let 
us now look at these equations in more detail when expressed in a global coordinate system (and not in 
local rest frame coordinates as done up to now). 


5 Transverse expansion dynamics in central collisions at high energy 

I will restrict my discussion here to azimuthally symmetric systems with longitudinal boost invariance. 
With this approximation we can describe the transverse expansion in central collisions at very high energy 
in a domain near midrapidity. A generalized discussion which relaxes the assumption of azimuthal 
symmetry and thus allows for collisions at any impact parameter has been recently given in 0 . 

Boost-invariant and azimuthally symmetric systems are conveniently described in (r, r, (j), rf) co¬ 
ordinates where r = \/t 2 —z 2 is longitudinal proper time, p = \ ln[(£+z)/(i— z)\ is space-time rapidity, 
and r = (r, 7) arc polar coordinates in the plane transverse to the beam direction z. Such systems are then 
characterized by macroscopic observables which are independent of o and //, and by azimuthally constant 
phase-space distributions which depend only on the difference Y—p (where Y = \ In [(E+p z ) / (E—p z )\ 
is the momentum-space rapidity of a particle with longitudinal momentum p z and energy E). We use 
lowercase latin letters to denote vector and tensor components in this curvilinear space-time coordinate 
system. The metric tensor in this coordinate system reads 

g mn = diag(l,-l,-l/r 2 ,-l/r 2 ), g mn = diag (1,-1, -r 2 , -r 2 ). (36) 


The flow velocity is parametrized as 


U m = 7 r (l, V r 


0 , 0 ) with y r = 


V 1 - v, 


(37) 


with radial transverse flow velocity v± = v r (r, r) e r and vanishing flow components lY and u n . For 
vectors and tensors, the usual Cartesian derivatives d ;t must be replaced by covariant derivatives, denoted 
by semicolons: 


d,,f -> j" = d m j 


+ C k j k 


f) rr it/A rpTll O rpTll , -pTT rpkl , rpnkj^l 

0^1 1 -m Um-L i 1 mk-L ' km’> 


(38) 


where T* fc = 
of r jk are 


7)fj' rn (Ojfjkrn — d m gjk) are the Christoffel symbols. The nonvanishing components 


= 

r]T rr) 


— v T = t = 4 -— 

T ’ VV ’ 4>r r<t> “ r > 


r 


= — r. 


(39) 
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The time derivative in the local comoving frame and the local expansion rate are thus computed as 

D = u ■ d = 7 r (ri T + v r d r ), (40) 

6 = d ■ u = — d T (Tj r ) + - d r ( Wr'yr ) (41) 

r r 

Due to azimuthal symmetry and longitudinal boost invariance, the n = (f> and n = p components 
of the equations of motion T mn - m = 0 are redundant. The n = r and n = r components can be written as 

El El 


p + n + T Z 7T 


\ B A rTTT )+\ d ’( rTT ") : 

-dr(rT Tr ) + -dr(r{T Tr v r + V r j) = + 


2 -ttVV 


p + n + r 2 TT^ 

r 

rjiTr rpT 


(42) 

(43) 


With the shorthand notations T mn = rTT mn , V r = rrV r , and v r = these are brought into 

“standard (Cartesian) form” 


d r T TT + d r (v r f TT ) = -r(p + n + r 2 7r^), 
d T T Tr + d r [y r T rr + Pr) = +T (p + n + r 2 7T^) . 

The coricsponding transport equations for the dissipative fluxes read 

[d T + v r d r ^j 
(d T + v r d r ^ i 

(d r + v r d r )u = —— [n + c^], 

with the following explicit expressions for the shear tensor components: 

■m 1 (0 7r' 


VV — 

1 

jrVV ^ 

'e 

lr\ 


7rT-7r 

71 r 2 

v3 

~ r)\ 


i 

Lh 2r ? 

(- 

'JrVr 



11 r 2 

\3 

r 


r z \S r J 


(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 


The hydrodynamic equations require the equation of state (EOS) p(e) for closure, i.e. after each 
transport step in time we must extract at each spatial grid point the boost velocity v r between the global 
and local rest frames and the local energy density e from the dynamical variables T TT and T Tr . The 
energy density is obtained from 

e = T TT -v r T rr , (51) 

where the radial velocity v r must be extracted from the implicit equation 


v r 


T 


rr 


T TT + p(e=T rT -v r T rr ) + n - r 2 ^ - t 2 tt w 


(52) 


by a one-dimensional zero search. This is still the same degree of numerical complexity as in the ideal 
fluid case f8); for dissipative hydrodynamics without azimuthal symmetry, however, this part of the 
problem becomes numerically more involved 0. 














6 First numerical results 


I close this talk by showing some preliminary results (7) from a numerical simulation of the equations 
derived in the preceding section, using a simple massless ideal gas EOS, p = |e, with e = aT 4 , a = (16 + 
T^/)fh' We ne §l ect bulk viscosity, £ = 0 . 

In classical kinetic theory, explicit expressions can be obtained for the viscosity coefficient // and 
the relaxation time r n in terms of the collision term. For a strongly coupled QGP, neither rj or tv are 
known. We treat them as phenomenological parameters. For guidance, we use perturbative lOMU and 
AdS/CFT-based ITU estimates for //, respectively, and a kinetic theory estimate |3) for tv. 

The shear viscosity coefficient ij for hot QCD was determined perturbatively to leading logarithmic 
accuracy in [10. dU. For a s ~ 0.5 the result in El gives ^ = 0.135. A lower limit for the shear 
viscosity in infinitely strongly coupled N = 4 SUSY YM theory and variations thereof was derived in 
GZ), exploiting the AdS/CFT correspondence: 2 = 0.08. In kinetic theory, in the Boltzmann gas 

approximation, the relaxation time is estimated as tv = 2 ///l 2 = 2 r/ -jb 0 

For the initial energy density distribution in the transverse plane, we used a Woods-Saxon param- 
eterisation, 

e(r) =-(53) 

1 + e o. 

with R = 6A fm, a = 0.54 fm. This is not very realistic, but facilitates comparison with the results of 
0 . (eo = aT 4 is the central energy density at initial time t = r,.) I show results for initial conditions 
Tj = 0.3 GeV and t, = 0.5 fm/c, with zero initial radial flow (v r (r, r,) = 0). 

For the non-ideal fluid, initial viscous pressures 7 r[T and are required. Even though v r and 
its derivatives are zero initially, due to the Bjorken longitudinal motion the stress tensor is not zero: 
T^a vv = — -y=r. r 2 (T 0ft> = We assume that at initial time r,;, the viscous pressure components are fully 
relaxed to the Bjorken scaling expansion values, 


rr 2 

7T- ■ = r 7T- ■ = 
ini ' mi 


T i ryn _ 

2 ini “ 3 ^ ' 


(54) 





Fig. 1: T/ = 160MeV surface for differ¬ 
ent initial viscous pressures. 


Fig. 2: Tf = 160 MeV surface for dif- Fig. 3: T/= 160 MeV surface for differ- 
ferent shear viscosities rj. ent relaxation times tv . 


Figures H to 01 show that the space-time evolution of non-ideal fluids depends sensitively on (i) the 
initial viscous pressure, (ii) the viscosity coefficient 77 , and (iii) the relaxation time T n . In all three figures 
we show the shape of the the Tf = 0.16 GeV surface, with the solid black line giving the ideal fluid result 
for reference. Viscous effects are seen to slow down the cooling of the matter, increasing the lifetime of 
the fireball and its average transverse size at freeze-out. 
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In Fig. |U the Tj = 0.16 GeV surface is shown for different initial viscous pressures, 7 r-T = 0, 
and respectively, using 2 = 0.135. The higher the initial viscous pressure, the more extended is 
the freeze-out surface and the larger are the deviations from ideal fluid dynamics. The life time of 
the dissipative QGP is extended by 20% if the initial viscous pressure is increased from zero to r^jr. 
The freeze-out surface also depends sensitively on the value of the viscosity coefficient (Fig. [2j>. As 
the viscosity decreases, the departure of the freeze-out surface from ideal behavior also decreases. In 
Flg-0we s how the freeze-out surface for different relaxation times, = 0.5 Tki n , ri, in and 2 r^;,, (where 
Tki n = 0), for fixed viscosity ri/s = 0.135. As the relaxation time is increased by a factor 4, the freeze- 
out time in the fireball center increases by 25%. 

All the viscous effects shown in Figures |j]to|3] increase if the initial time t, is decreased, keeping 
the total fireball entropy and all other parameters constant. This is due to the increasing initial longi¬ 
tudinal expansion rate 6 = — which results in a larger ratio —. This ratio is the figure of merit which 
controls the importance of viscous corrections to ideal fluid dynamics. Our studies show that there is 
not only a minimum thermalization time r t h erm ni ter which ideal fluid dynamic can be applied, but there 
is also a minimum time t, < r t h e rm for the applicability of viscous fluid dynamics. The initial condi¬ 
tions for viscous hydrodynamics at that time r t must be obtained by matching the decomposition ifTTt of 
the energy-momentum tensor to the corresponding result for T IW from some preceding non-equilibrium 
kinetic evolution. 
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